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ABSTRACT 

We report the detection of thermal emission from the hot Jupiter WASP-3b in 
the Kg band, using a newly developed guiding scheme for the WIRC instrument 
at the Palomar Hale 200in telescope. Our new guiding scheme has improved 
the telescope guiding precision by a factor of ~5-7, significantly reducing the 
correlated systematics in the measured light curves. This results in the detec- 
tion of a secondary eclipse with depth of 0.181%±0.020% (9-er) - a significant 
improvement in WIRC's photometric precision and a demonstration of the capa- 
bility of Palomar/WIRC to produce high quality measurements of exoplanetary 
atmospheres. Our measured eclipse depth cannot be explained by model at- 
mospheres with heat redistribution but favor a pure radiative equilibrium case 
with no redistribution across the surface of the planet. Our measurement also 
gives an eclipse phase center of 0.5045±0.0020, corresponding to an ecoso; of 
0.0070±0.0032. This result is consistent with a circular orbit, although it also 
suggests the planet's orbit might be slightly eccentric. The possible non-zero ec- 
centricity provides insight into the tidal circularization process of the star-planet 
system, but also might have been caused by a second low-mass planet in the sys- 
tem, as suggested by a previous transit timing variation study. More secondary 
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eclipse observations, especially at multiple wavelengths, are necessary to deter- 
mine the temperature-pressure profile of the planetary atmosphere and shed light 
on its orbital eccentricity. 

Subject headings: Infrared: planetary systems - Planetary systems - Stars: 
individual (WASP-3), 
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from ground for about a dozen hot Jupiters since 2009 (e.g.. 



Alonso et alJ20Q9: 



2011 



Zhao et al. 



2012 



Gibson et alJBOld : iCroll et al.ll2010aL 12011a ; 
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Sine & LoDez-Morales! 


2009c 


Caceres et al. 


2011; 


de Mooii et al. 



etc.). These observations provide important probes of planetary at- 
mospheres at the near-IR where most of the bolometric output of the planet emerges. They 
also measure deeper and higher-pressure layers of atmospheres than observations at longer 
wavelengths, highly complementary to the Spitzer IRAC measurements. When combined 
together, these broad-band multi-wavelength measurements can provide constraints to plan- 
etary SEDs, and help to dis tinguish between differing atmospheric pressure-temperature 



profiles and chemistries (e.g., iMadhusudhan et al.ll201ll ). Detection of secondary eclipses of 



transiting planets also provide important constraints to their small orbital eccentricities that 
are usually hard to be distinguished from zero by radial velocity me asurements. Precise ec- 
centr icity measurements can improve estimates of planetary radius (IMadhusudhan fc Winn 
20091 ). and can provide important information to their tidal circularization process, shed- 
ding light on the nature of, in some cases, anomalously inflated radii of some planets (e.g., 
Miller et alJboogh . 



WASP -3b is a massive tran siting hot Jupiter (1.76Mj up ) orbiting a F7-8V type star at 
0.0317AU ( IPollacco et al.l 120081 ). Several groups have measured its spin-orbit alignment via 
the Rossiter-McLaughlin effec t, finding a close alignment between the stellar rotation axi s 
and the planet's orbital axis (ITripathi et al.l l2010t iMiller et al.l 12010c ISimpson et al.l 120101 ) . 
Studies have also measured transits of WASP-3b at many epochs to sea rch for possibl e tran- 
sit timing variation (TTV) caused by a low-mass body in an outer orbit. iGibson et al.l (120081 ) 
observed two transits of WAS P -3b but did not find sign ificant difference from the original 
timing of IPollacco et al.l (120081 ). iMaciejewski et al.l (120101 ) later combined 6 new transits with 
previous data and found a periodic time variation of ~ 3.7 days, which could be interpreted 
as caused by a hypothetical second planet with a mass of ~15 M e at a semimajor axis of 
0.0507AU and a period very close to the 2:1 mean motion resonance. However, they also 
emphasized that more observations are required to confirm this periodic variation. Mean- 
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while, iLittlefieldl (1201 If) observed 5 tr ansits of WASP-3b and found supportive evidence to 
the result of iMaciejewski et al.l ( 120101 ) . Efforts have also been employed to se arch for addi- 
tiona l transiting planets in the WASP-3 system, but no candidates were found ( [Ballard et al. 
201lh . 



The close-in orbit and relatively large radius of WASP-3b and the strong radiation from 
its host star (T e ff ~6400K) make its atmosphere highly irradiated (T e(? =1960K), making 
it an ideal target for secondary eclipse detections and studies of its atmospheric properties. 
Despite its high temperature, WASP-3b's atmosphere remains one of the least characterized 
among the most irradiated hot Jup iters, i.e., with only an upper limit o f secondary eclipse 



at 650nm (IChristiansen et al.l 1201 ll ) and two detections in the Ks band ( jCroll et al.l 1201 lb 
Croll et al. 2012 in preparation). 

Here we report another detection of WASP-3b's thermal emission in the Ks band us- 
ing the Palomar Hale 200in telescope with an improved guiding scheme. We present our 
observations, including the guiding improvement, and data reduction procedure in §2j We 
describe our data analysis and results in §3J We discuss our measured eclipse phase center 
and compare the eclipse depth of WASP-3b with existing models in §4J We then summarize 
our results in SJ5l 



2. Observations and data reduction 

The observation of WASP-3b was conducted in the Ks band with the WIRC instrument 



(IWilson et al.l 120031 ) on Palomar 200-in Hale telescope on UT 2011 June 24 (PI: Hinkley). 
WIRC has a 2048 x 2048 Hawaii-II HgCdTe detector with a scale of 0.2487"/pixel and a 
wide field of view of 8.7' x 8.7'. The observation started at 04:05:59.012 UTC on 2011 June 
24, and ended 358.12 minutes later. To minimize instrument systematics, we stayed on the 
target without dithering for the entire observation. The telescope was defocused to a FWHM 
of about 2.5" - 3" to keep the counts well within the linearity regime and to mitigate intra- 
pixel variations. Each image was taken with 12sec exposure and single fowler sampling. A 
total of 683 images were obtained. The duty cycle of the observation was 44%. 



2.1. Improved guiding for Palomar /WIRC 

Because WIRC does not have a dedicated guider, its guiding of targets relies completely 
on telescope tracking. A previous study has shown that the limited tracking precision of the 
telescope, especially along the X-axis (i.e., the direction of R.A.), results in highly correlated 
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systematics in the detected light curves (jZhao et al.ll2012l ). These systematics are caused by 



inter-pixel variations of the detector (due to imperfect flat fielding) and is a dominant source 
of "red noise" in high precision light curve measurements. Stabilized guiding is necessary to 
partially mitigate this problem. We therefore designed an active guiding scheme to correct 
for the telescope tracking errors during observing based on the information obtained from 
previous images. The algorithm has been integrated into the WIRC control system to offer 
fast corrections. Figure [1] demonstrates the effectiveness of the new guiding scheme, which 
has improved the guiding precision of WIRC by a factor of ~5-7. Currently, the precision is 
limited to ~2-3 pixels due to a highly periodic gear oscillation with a frequency of ~0.5Hz 
from the telescope. We have been developing another algorithm to further correct this 
oscillation, and expect to deliver better guiding precision in future observations. 

The observation of WASP-3 was conducted during the development stage of the new 
guiding scheme. The guiding performance was similar to the middle panel of Figure HJ 
despite a loss of 45.43min of data due to a software glitch occurred during middle eclipse 
(also see Figure [2]) and a ~ 0.5-pixel centroid shift after the observation was recovered. 



2.2. Data reduction 

For the data reduction process, we first subtracted all images with corresponding av- 
eraged dark frames. Twilight and sky flats were normalized and averaged to get a master 
flat field. A bad pixel mask was then created with the master flat and dark frames. The 
bad pixels in each image were interpolated with cubic splines based on adjacent flat-fielded 
pixels. WASP-3 is the brightest star in its relatively sparse field. To properly correct the 
highly correlated common-mode systematics in its light curve, six well separated and evenly 
distributed stars within the flux range of 0.13 to 0.71 times of that of WASP-3 were selected 
as references. Other stars in the field were too faint to have sufficient signal-to-noise and 
were thus excluded. We calculated the centroids of all stars in each image using a center- 
of-mass calculation, since it provided the smallest scatters of their relative positions. The 
time series of WASP-3's centroid was determined by averaging the relative positions of all 
reference stars after correcting for their relative distances. 

Aperture photometry was performed on WASP-3 and the reference stars following the 
IDL routines of DAOPHOT. The extracted fluxes of each star were normalized to the median 
of the time series. We used the median of the 6 reference time series as the final reference 
light curve due to the presence of outliers in the light curves of some reference stars and 
the fact that median is a more robust estimator. The final reference light curve is then 
used to normalize the flux of WASP-3 to correct for the common-mode systematics such as 
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variations of atmospheric transmission, change of seeing and airmass, etc. We applied 48 
different aperture sizes with a step of 0.5 pixel, and determined that an aperture with a 
radius of 16 pixels (i.e., a 32-pixel diameter) gives the smallest out-of-eclipse and in-eclipse 
scatters for the normalized WASP-3 data; this was taken as the final photometry aperture 
for all stars in every image. A sky annulus with 30-pixel inner radius and 35-pixel width was 
used for background estimation. The median value of the sky annulus was then used as the 
final sky background for subtraction. We have also explored different annulus ranges and 
sizes, and found consistent results. The top two panels of Figure [2] show the reduced fluxes 
of all 7 stars and the final normalized flux of WASP-3, respectively. The UTC mid-exposure 
time of each image w as converted to BJD^db using the UTC2BJD routine provided by 



Eastman et al.l (12010). T he orbital phases were calculated based on the latest ephemeris of 



Christiansen et al.l ( 120111 ). since they have taken both the previously published transit times 
and their new EPOCh transits into account (i.e., period=1.8468373±0. 0000014 days, and 
transit epoch T ( J BJD TDB )=2454686.82069±0.00039). 



3. Analysis and results 



After normalizing the time series of WASP-3 with the reference light curve, the flux 
decrease caused by WASP-3b's eclipse becomes visually identifiable (Figure To measure 
the eclipse depth and determine the phase center, we fit a light curve simultaneously with 
a background baseline to the data. Thanks to the improved guiding of WIRC, the drift of 
centroid is relatively small and we do not see obvious correlation between flux and centroid 
positions (Figure [3]). We experimented with both a linear baseline and a quadratic baseline 
with the data, and a linear baseline with a neglig ible slope is p referred to a quadratic baseline 
by the Bayesian Information Criterion (BIC) □ ( Liddld 120071 ). 



The light curve is generated following the prescription of iMandel fc Agoll ( 120021 ). assum- 
ing uniform bodies without limb-darkening. The stellar and planetary parameters for the 
l ight c urve (R p , R s tar, inclination, and semimajor axis) are adopted from IChristiansen et al. 
( 1201 it ). The free parameters in the least-square fit are: the eclipse depth, the mid-eclipse 
phase, the level of the out-of-eclipse baseline a±, and the baseline slope a<i- The known 
durations of ingress and egress are maintained in the fit. 



We employed the Levenberg-Marquardt (LM) algorithm (jPress et al.lll992l ) for the least 



1 We use a linear baseline of the form: / = a\ + a% • t, where / is the flux, t is the time of each datum, 
and m, ct2 are the linear coefficients of the baseline. The linear baseline gives a BIC value of 710, less than 
the value of 716 from the quadratic baseline. Thus, the linear baseline model is preferred. 
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square fit. To ensure that we find the global minimum instead of local minima, we searched 
the parameter space extensively with a fine grid of starting points on top of the least-square 
fit. The grid has a few hundred steps for each parameter. The fact that most starting 
values on the grid converge to the same minimum suggests that we indeed have found the 
global minimum. The data points are uniformly weighted such that the xl is nearly 1.0. 
The global best-fit light curve gives an eclipse depth of 0.181%±0.020%, and a phase center 
of 0.5045±0.0014. The best-fit model is shown in Figure [2j along with the residual of the 
best-fit and the binned light curve. The right panel of Figure [3] compares the noise level 
of WASP-3b's light curve with the Gaussian noise expectation. Both the in-eclipse and 
out-of-eclipse data follow the Gaussian expectation closely, although there are still some 
uncorrected systematics in the data. 

To examine the statistical significance and robustness of the eclipse depth and to es- 
timate its error, we conduct 2 statistical tests. We first apply the standard bootstrapping 



technique ( jPress et al.lll992l ). In each bootstrapping iteration, we uniformly resample the 
data with replacement. For each new sample, we re-fit the light curve and baseline model to 
determine the eclipse depth and phase center, using the aforementioned grid search and LM 
minimization. This technique is suitable for unknown distributions like our case, and can 
robustly test the best-fit model and the distribution of the parameters. A total number of 
2000 iterations are performed and the resulting distributions of the eclipse depth and phase 
center are nearly Gaussian, with a median depth of 0.183%±0.019% and an eclipse phase 
center of 0.5045±0.0020, highly consistent with the previous best-fit. 



For the second test, we use the "prayer-bead" residual permutation method (IWinn et al. 



20081 . and references therein). The same light curve and baseline model are employed to 
re-fit the simulated data in each iteration. A total number of 1365 iterations (i.e., 2N-1, 
where N=683 is the number of data points) are conducted. This method maintains the 
time-correlated errors and is therefore another robust way of testing our fit. Thanks to 
the minimal amount of "red noise" in the data, the resulting distribution is also close to 
Gaussian. The resulting median depth and 1-er error is 0.185%±0.019%, and the eclipse 
phase center is at 0.5044±0.0015, also consistent with the previous results. 

We take the values from the LM best-fit and the largest error bars from the 3 different 
error analyses as our final results, and summarize them in Table [TJ Based on the average 
flux of the target and the sky background, the expected photon noise precision for the eclipse 
data is 0.0065%. Thus, our precision of 0.020% corresponds to ~3 times of the photon noise 
limit. 
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Discussion 



We compare our K q band measur e ment of W ASP- 3b in Figur e B] with atmospheric mod- 



els generated base d on 



Barman et al. 



(120051 ) and iBarmanl (120081 ). Our measured flux ratio 
agrees with that of ICroll et all (l2011bl ). 0.176%l{j;{$?. However, it is too high to match mod- 
els with heat redistribution. Instead, the flux ratio is more consistent with a hot atmosphere 
in pure radiative equilibrium (at the 2-cr level), with no redistribution across the surface 
of the planet. Such a model has a nearly isothermal radial temperature profile across the 
near-IR photosphere over most of the dayside, except near t he planet limb. Th e nightside 
is very cold. This conclusion is also consistent with that of ICroll et al.l (l2011bl ). Perhaps 
coincidentally, the observed flux matches very closely the value predicted by a model for 
the substellar point (top curve), suggesting that, at the K^-band photospheric depth, the 
temperature is on average close to that at the substellar point (also corresponding to a black 
body temperature of ~2435K). It is unrealistic for a planet's entire dayside to be as hot as 
the substellar point, indicating a large departure from radial temperature profile predicted 
by traditional one- dimensional atmosphere models. However, with no color information, it 
is difficult to infer detailed information about the nature of the temperature structure. More 
observations at other wavelengths are definitely needed. 



Our final measurement of the eclipse center, 



0.5045 ± 0.0020 corresponds to a 



delay of 11.97±5.32 minutes from the expected mid-eclipse time based on lChristiansen et al. 



2011 



(j2008|) 



Pollacco et al. 



This gives ecosa; = 0.0 070 ± 0.0032, c onsist ent with the values of 
e = 0.05 ± 0.05), and Isimpson et all (boioh ( e = 0.07 ± 0.08). However, this 
result is only consistent with the joint measurement of ICroll et al.l ( 1201 lbl ) at 2.5-cr level 
(cj) = 0.4999tQ;Qoio) ; although it agrees better with their sec ond and higher-sign i ficanc e mea- 



surement of <j> = 0.5014+°;°°^. In search for TTV signals, iMaciejewski et al.l ( 120101 ) founc 



a potential periodic time variation with a semi-amplitude of 0.0014 days, which could be 
interpreted as a hypothetical low-mass body in an outer orbit. They also made a joint 
reanalysis of the existing RV data and found an eccentricity of 0.05 ± 0.04, which is also 
consistent with our slight non-zero result. However, they also pointed out that the possible 
non-zero ec centricity might be a result of confusion with a two-planet system in an inner 2:1 



resonance. 



Miller et all (120091 ) indicated that the observed radius of WASP-3b, 1.385R Jup , 
is more consistent with their full tidal evolution model. Thus, if the orbit of WASP-3b is 
indeed slightly eccentric, it might be able to shed light on its migration and tidal evolution 
history. However, we also emphasize th at the current mea surement is still consistent with 
zero and given the slight difference from ICroll et al.l (j2011bl ) , more observations are required 
to better constrain the eccentricity and shed light on the aforementioned scenarios. 
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Conclusions 



We have developed a new integrated guiding scheme for Palomar/WIRC observations. 
Our algorithm has improved the guiding precision by a factor of ~5-7, significantly miti- 
gated the centroid drifts of the targets on the detector, and largely reduced the correlated 
systematics in the light curves seen in our previous study Using the new guiding scheme, 
we have detected the K$ band thermal emission of the hot Jupiter WASP-3b at 9-a signif- 
icance. The detected secondary e clipse h as a depth of 0.181%±0.020%, in agreement with 
the previous result of ICroll et al.l ( 1201 lbl ). The measured flux ratio of the planet is too high 
to be explained by models with heat redistribution but favors a pure radiative equilibrium 
case with a very cold nightside. Further observations at multiple wavelengths are necessary 
to help determine the temperature-pressure profile of the planetary atmosphere and shed 
light on the nature of its high eclipse depth in Kg- Our measurement also gives an ecosw 
of 0.0070±0.0032, consistent with a circular orbit while also sugge sting; the pl a net's o rbit 
might be slightly eccentric. This result differs sightly from that of ICroll et al.l ( l2011af ) by 
~2.5-cr. On the other hand, a previous study has found possible periodic TTV signals in 
the syst em, and a small non-ze ro eccentricity might be caused by a second planet in the 
system ( jMaciejewski et al.ll2010l ). More secondary eclipse observations are certainly needed 
to better constrain the eccentricity. 

Despite the 45-minute gap during the eclipse, we still achieved a 9-a detection thanks to 
the large aperture of the telescope and the substantially reduced systematics. This demon- 
strates the capability of Palomar/WIRC in providing high quality secondary eclipse mea- 
surements of hot Jupiters, and its potential to expand to other wavelengths such as H and 
J. With the large aperture size and additional improvements in guiding, we expect Palo- 
mar/WIRC to make significant contributions to the studies of exoplanetary atmospheres. 
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Fig. 1. — Comparison of the centroid drifts in the x direction (R.A.) before and after im- 
plementing the new guiding scheme. The top panel shows the centroid drift with a speed of 
about 5-10 pixels per 20min (0.75-2.5arcsec), resulting in large systematics in the data of 
2011Jan29 (IZhao et al.ll2012l ). Abrupt jumps of the centroid during that night were due to 
manual corrections of telescope position. After applying the guiding algorithm to WIRC on 
UT 2011Jul27 (middle panel), the telescope drift was mostly corrected, leaving a residual 
fluctuation of only ~2 pixels stemming from a telescope gear oscillation. The observation of 
WASP-3 on UT 2011 June 24 was carried out during the development stage of the scheme 
and is shown in the bottom panel. 



-12- 



1.02 



0.98 

0.96 
0.94 



£ 1.005 

13 
CD 

1.000 

nj 
E 

o 0.995 



0.010 



0.005 




CD 

en 



0.000 
-0.005 - 

-0.010 
1.002 

x 1.001 

f 1.000 

g, 0.999 

| 0.998 

< 0.997 

0.996 







Best-fit Depth = 0.181 +/-0.020% 



Nsum=18 



46 0.48 0.50 0.52 0.54 0.56 0.58 0.60 
Phase 



Fig. 2. — Reduced flux and best-fit light curve of WASP-3b. The first panel shows the 
reduced and normalized flux of WASP-3b (black dotted line), overplotted with the fluxes of 
the 6 reference stars (colored lines). The second panel shows the light curve of WASP-3b 
after correcting with the reference light curve, along with its best-fit model (solid blue line). 
The third panel shows the residual of the best-fit. The bottom panel shows the averaged 
light curve together with the best-fit model. Error bars of the points are calculated from the 
scatter of the data in each bin. A software glitch occurred during middle transit, causing a 
gap of 45.43min in the light curve. The larger scatter in the data between phase 0.54 and 
0.55 was likely due to deteriorated seeing. 
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Fig. 3. — Left: Flux of WASP-3b as a function of x and y positions of its centroid. The 
centroid drifts in both directions are relatively small and have no obvious correlations with 
flux. Right : Comparison of WASP-3b's noise level with Gaussian noise expectation. Both 
the in-eclipse and out-of-eclipse data follow the Gaussian expectations closely as the data 
points are binned down. 



Table 1. WASP-3 eclipse parameters 



Parameter Final result 



& 0.181%±0.020% 

(pmid-eclipse 0.5045±0.0020 

ecosw 0.0070±0.0032 



Adopted parameters 1 



T 2454686.82069 (BJD) 

P 1.8468373 (days) 

R P 1.385 (Rj U p) 

R» 1.354 (R Q ) 

a 0.03167 (AU) 

i 84.22° 



^om lChristiansen etaD l|201ll l 
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Fig. 4. — Comparison of Ks measurement with atmospheric models. The blue dotted line 
shows the planet- wide full redistribution model (also known as the 4n redistribution). The 
green line shows the dayside-only redistribution model (2n redistribution), the orange line 
shows the model with no heat distribution at all, while the red line shows the flux emerging 
from the planet assuming everywhere is as hot as the substellar point. All models have an 
inverted temperature profile for P gas <1 to 10 bar. The atmosphere becomes isothermal 
at larger P gas . The data point from this work is shown in black dot, along with the filter 
transmission profile (black line at the bottom). The band-averaged model points are shown 
as colored diamonds. 



